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Abstract 

The measurement of muon energy is critical for many analyses in large Cherenkov detectors, particularly those that 
involve separating extraterrestrial neutrinos from the atmospheric neutrino background. Muon energy has traditionally 
been determined by measuring the specific energy loss (dE/dx) along the muon's path and relating the dE/dx to the 
muon energy. Because high-energy muons (E^ > 1 TeV) lose energy randomly, the spread in dE/dx values is quite 
large, leading to a typical energy resolution of 0.29 in log 10 (£' iU ) for a muon observed over a 1 km path length in the 
IceCube detector. In this paper, we present an improved method that uses a truncated mean and other techniques to 
determine the muon energy. The muon track is divided into separate segments with individual dE/dx values. The 
elimination of segments with the highest dE/dx results in an overall dE/dx that is more closely correlated to the 
muon energy. This method results in an energy resolution of 0.22 in log 10 (£ , /i ), which gives a 26% improvement. This 
technique is applicable to any large water or ice detector and potentially to large scintillator or liquid argon detectors. 

Keywords: 

muon energy, dE/dx, neutrino energy, truncated mean, Cherenkov, IceCube detector 
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1. Introduction 

Large ice or water Cherenkov detectors may observe 
up to 100,000 Vp, per year yj]. These events are used 
for a wide variety of analyses, including searches for 
point sources of neutrinos J2, S 01, diffuse extrater- 
restrial neutrinos fl. S Si, standard and non-standard 
neutrino oscillations [7], and measurements of the total 
neutrino-nucleon cross-section via neutrino absorption 
in the Earth. 

These analyses rely upon the measurement of the 
neutrino energy, which is determined from the energy 
of the muon that is created in the neutrino interaction. 
Above ~ 1 TeV, the muon energy is usually determined 
by measuring the specific energy loss, dE/dx, of the 
muon as it travels through the detector. The Cherenkov 
photons from the muon and also those derived from 
the charged particles produced by stochastic (random) 
muon interactions are then detected. This approach is 
disadvantageous because, for > 1 TeV, muons lose 
most of their energy stochastically, and a small number 
of high-energy interactions will not only skew the mean 
but also enlarge the spread in dE/dx values. 

In this paper, we present an improved method for 
calculating the muon energy loss, which leads to sig- 
nificant improvement in the energy resolution. Instead 
of averaging the muon dE/dx over the entire observed 
muon path length, we divide the path into independent 
segments, or bins. The dE/dx is calculated separately 
for each bin. Then, the bins with the highest dE/dx 
values are discarded before calculating a new average 
dE/dx, thus producing a truncated mean. This method 
is successful because the truncated mean minimizes the 
effects of the large stochastic events which would other- 
wise skew the mean and enlarge the spread. 

The truncated mean method has previously been cal- 
culated for muons fl, although that analysis did not 
use it to determine the muon energy. The method was 
also explored at a basic level for the DUMAND project 
J^tl . The method has parallels to the one that was used 
to identify pions, kaons, and protons in wire chambers 
from their specific energy loss dE/dx in the gas. By dis- 
carding the highest 30% or 40% of the dE/dx measure- 
ments from the wires, the energy resolution was greatly 
improved 1 10, 11]. As with wire chamber dE/dx mea- 
surements, the muon energy resolution is improved by 
discarding the most energetic stochastic losses. This 
method should be applicable to any large water or ice 
Cherenkov detector, such as IceCube Oil, ANTARES 
IH, or the proposed KJVBNeT |l3 
or Hyper-Kamiokande (Hyper-K) 



gon detectors, such as the Low Energy Neutrino Astron- 
omy detector (LENA) |16] or the Long-Baseline Neu- 
trino Experiment (LBNE) ll7Tl. 



2. Muon Energy Loss 

Muons lose energy via ionization and by stochastic 
processes such as bremsstrahlung, pair production, and 
photonuclear interactions H18L I19L l2Qfl . Muons also emit 
Cherenkov radiation, but this is a very small fraction 
of the total energy loss. Ionization loss is roughly in- 
dependent of muon energy and essentially constant per 
unit length, while the dE/dx due to stochastic processes 
rises linearly with the energy. Thus the total average 
energy loss of the muon is 112 If 



dEf, 
dx 



= A + BE,, 



(D 



MEMPHYS 114(1 , 
1511 . This approach 
may also be useful for proposed scintillator or liquid ar- 



where A as 0.0024 GeV (g/cm 2 ) -1 accounts for the en- 
ergy loss due to ionization, and B w 0.000032 (g/cm 2 ) -1 
is due to the stochastic energy loss l20ll . 

Two main methods are used to measure E^, For 
lower-energy muons, the energy is proportional to the 
path length. This approach works if the muon is con- 
tained within the detector, and the starting and stopping 
points can be determined. For example, a 300 GeV 
muon has a most probable path length of ~1 km in ice. 
The most probable path length is ~2 km for a 1 TeV 
muon, and the length rises only logarithmically with E^, 
reaching 20 km for a 1 PeV muon ll20ll . 

At energies above a few hundred GeV, one measures 
the dE/dx and infers E M from that value, either by us- 
ing Eq.Q]as an approximation or by calculating the cor- 
relation using Monte Carlo simulation. The total en- 
ergy loss in the entire detector is summed, either us- 
ing calorimetric measurements of dE/dx (such as by 
the Frejus collaboration) i22ll or by observing the light 
produced by the stochastic interactions of the muon and 
charged particles lEill . Photonuclear interactions pro- 
duce slightly less light than electromagnetic interac- 
tions such as bremsstrahlung and pair production, but 
the photonuclear fraction is well known, so this can be 
compensated for. 

The IceCube Neutrino Detector is highly suited for a 
study involving muon energy loss. This analysis used 
IceCube in its 2010 configuration, when the detector 
consisted of 79 vertical strings (with roughly 125 m 
horizontal spacing in a triangular grid), with 60 optical 
modules per string (with 17 m vertical spacing) between 
1450 and 2450 m below the surface of the ice sheet at 
the South Pole. This included DeepCore, which is a 
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denser subarray of strings near the detector center with 
more-sensitive optical modules placed preferentially on 
the bottom 400 m of each string 12411 . Each optical mod- 
ule holds a 25 cm photomultiplier tube, plus associated 
electronics for data acquisition, control, and calibration 

a. 

Simulated muon trajectories were determined using 
a program called neutrino-generator (NUGEN) which 
is based upon the All Neutrino Interaction Simulation 
(ANIS) ||25l]. NUGEN simulation includes several ef- 
fects, including the ice/rock boundary below the detec- 
tor, Earth-neutrino absorption using the density profile 
from the Preliminary Reference Earth Model (PREM) 
i26ll . neutral current regeneration, etc., but it does not 
include neutrino oscillations. The program starts with 
an isotropic distribution of neutrinos on the Earth's 
surface, following, in this case, a E~ l spectrum for 
10 GeV < E Y < 1 EeV. These neutrino events can be 
reweighted to simulate softer spectra. NUGEN then 
propagates the neutrinos to the South Pole, account- 
ing for neutrino absorption in the Earth. Finally, NU- 
GEN simulates the neutrino interactions. The result- 
ing muons are propagated through the ice using the 
Muon Monte Carlo (MMC) code [20], which includes 
detailed models for muon energy loss and the Landau- 
Pomeranchuk-Migdal effect 127 1. 

From the energy deposition determined by MMC, a 
program called Photonics determines the number and 
timing distribution of expected photoelectrons for a 
given optical module 12811 . Photonics also accounts 
for the detector geometry of the optical modules. Fi- 
nally, a set of detector simulation programs models the 
response of the photomultiplier tube, electronics, and 
triggers lf29tl . Then the simulated data is analyzed us- 
ing standard IceCube reconstruction programs. Muon 
tracks are typically reconstructed with an angular reso- 
lution of < 1 degree ll30Tl and a positional accuracy (in 
the plane perpendicular to the track) of a few meters. 

2.1. Conventional dE/dx 

Before discussing the truncated mean dE/dx in de- 
tail, we will discuss the more conventional approach. 
First, the number of photoelectrons observed by all of 
the optical modules is summed. Then, the expected 
number of photoelectrons for a fixed energy loss of 
1 GeV/m with the same trajectory is calculated. Since 
the photoelectron yield is expected to be directly pro- 
portional to the muon energy loss rate, the calculated 
dE/dx value is equal to the ratio of the observed num- 
ber of photoelectrons to the expected number for the 
1 GeV/m track, multiplied by 1 GeV/m. This linear rela- 



tionship is a proper dE/dx approximation for muon en- 
ergies above ~ 1 TeV. Because the method uses the dis- 
tance between each sensor and the muon track, it auto- 
matically accounts for varying track geometries, sensor 
densities, and optical properties of the detector medium, 
since those effects are included in both the observed and 
expected numbers of photoelectrons. 

Figure Q] shows a scatter plot of simulated muon en- 
ergy versus calculated dE/dx. Because the entire track 
is used, the simulated muon energy value is extracted 
from the track's closest approach to the center of the de- 
tector (instead of at the entrance to or exit from the de- 
tector volume). The dE/dx values have a large spread. 
Some of the spread is due to the limited sampling by the 
sensors and consequently the uncertain measurement of 
the energy loss. However, most of the spread is from the 
event-by-event variation in stochastic energy losses. 



10 E 




12 3 

log (Untruncated dE/dx [GeV/m]) 



Figure 1: (color online) Simulated vs. conventional (untruncated) 
dE/dx at the center of the detector. The dE/dx value is calculated 
from the ratio of the observed energy deposition to the expected de- 
position for a minimum ionizing muon with the same trajectory. The 
black solid line uses Eq.[T]for the "linear" fit. The blue dashed line 
uses the 3-equation fit described by Eq.[2]and Eq.f5] 



The spread in dE/dx values for a given true energy is 
quite large for the conventional (untruncated) method. 
Figure [2] shows a typical example of the spread for a 
narrow energy slice, with this one at ~10 TeV. The ex- 
cess of events to the right of the Gaussian curve is a 
direct result of the large stochastic losses. 

The desired output of the method is the energy for 
each muon. This is done using a fit to the curve in Fig.Q] 
to calculate from the dE/dx value, and then we can 
compare the calculated energy to the simulated energy. 

It is apparent that the "linear" fit curve from Eq. Q] 
does not follow the data precisely, particularly at higher 
energies. The main experimental effect is saturation in 
the photomultiplier tubes for very high light intensities 
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2 3 
log (dE/dx [GeV/m]) 



Figure 2: (color online) Conventional (untruncated) muon dE/dx, for 
muons with simulated energies ~10 TeV (logjofE^) = 4.0 ± 0.05). 
The black line is a Gaussian fit, with cr = 0.29 ± 0.01, and RMS 
= 0.44. The excess of events to the right of the Gaussian curve are 
muons with large stochastic losses. 



03111 . This saturation was included in the simulations, 
but it is not corrected for in the reconstruction software. 
This produces the upturn in the curve at higher ener- 
gies for dE/dx > 10 3 GeV/m. For dE/dx < 1 GeV/m, 
there is also a small downturn in the curve, making it 
S-shaped. This occurs because the muon range at lower 
energy is less than the size of the detector, while the 
truncated method assumes that the track is infinite. 

We examined a variety of simple fit equations to im- 
prove our ability to characterize the muon energy res- 
olution. The optimal fit had three separate sections: 
a second-order logarithmic (log 10 ) polynomial fitted 
to the lower energy curve (for dE/dx < 1.5 GeV/m) 
and another second-order logarithmic polynomial to the 
higher energy curve (dE/dx > 500 GeV/m): 



log^^A' + B'loggj + C'log 2 ^ 



(2) 



where E M is the energy of the muon (in GeV) at the 
track's closest approach to the center of the detector, 
dE/dx is the calculated value (in GeV/m), and A', B', 
and C are dimensionless constants from the best fits to 
the curves. For the mid-range, we used the logarithmic 
version of the linear dE/dx equation (Eq.[T]i: 



log<£„> = log l^; 



dE A' 
dx 



(3) 



where is the energy of the muon (in GeV), dE/dx 
is the calculated value (in GeV/m), and A" and B" are 
dimensionless constants from the best fit to the curve. 
Because of PMT saturation and other detector effects, 
the constants in Eq.[3]must be adjusted for a proper fit. 



The new fit curve is shown in Fig. [T] as a blue dashed 
line. We used this 3-equation fit throughout this analy- 
sis, instead of the linear fit. 

The energy residual E res is nominally defined as 
E reco - E S j m , where E reco and E S i m are the recon- 
structed and actual simulated energies, respectively, 
at the track's closest approach to the center of the 
detector. Since this analysis covered a wide range 
of energies, the energy residual was redefined as 
E res = log 10 (£ reco ) - log m (E sim ). Figure [3] shows the 
results of this energy comparison, using the 3-equation 
fit. The distribution is skewed, with a long tail of events 
containing large stochastic energy deposition, indicat- 
ing that the energies for these events are badly overesti- 
mated. About 5.4% of the muon energies are overesti- 
mated by more than a factor of 5. 



a-2000 



O1500 
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Figure 3: (color online) Conventional (untruncated) muon energy res- 
olution, using the 3-equation fit curve from Fig. [T] for muons with 
simulated energies between 1 TeV and 1 EeV. The black line is a 
Gaussian fit, with the peak offset by -0.08 in log[o(£,,), cr = 0.29, 
and with an RMS of 0.37. The excess of events to the right of the 
Gaussian curve is a direct result of the large stochastic losses. 

If the distribution of energies is fit to a Gaussian, the 
energy resolution cr is 0.29 in log I0 (£ , /J ), or roughly a 
factor of 2.0 in energy. The statistical error is less than 
1 % for this resolution and all other resolutions presented 
in this paper (unless otherwise noted). However, the 
peak of the curve is offset from by -0.08, and the curve 
is skewed. The shift towards lower muon energies oc- 
curs because the fit curve in Fig.[T]does not line up with 
the highest concentration of events. This occurs because 
the fit curve is pulled towards higher dE/dx values for 
a given muon energy due to the stochastic losses. For 
a specific dE/dx value, the calculated energy using the 
fit equations for the bulk of the events will be too low. 
This pushes the peak of the distribution of energy resid- 
uals toward the negative. Fortunately, in a full analy- 
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sis, this bias would cancel out if the method were used 
identically for both data and simulation, and the spread 
would be the most relevant issue. 

In this analysis, we used both the Gaussian cr and the 
RMS to characterize the improvements in the spread in 
calculated dE/dx and E^, In most cases, the RMS was 
significantly larger than cr because the RMS is more 
sensitive to the high-energy tails in the distributions. 

2.2. Truncated Mean dE I dx 



The truncated mean method divides the track into 
bins which are bordered by planes perpendicular to the 
track, and the sensor data are binned accordingly by 
their location. The dE/dx value is determined as be- 
fore by finding the ratio of the observed photoelectrons 
to the expected photoelectrons, but a separate ratio is 
determined for each bin instead of the event as a whole. 
Then a fraction of the bins (using the optimization dis- 
cussed below) with the highest ratios are discarded, and 
the truncated dE/dx is calculated by summing the re- 
maining observed photoelectrons and expected photo- 
electrons and creating a new ratio. 

Many factors influence the choice of bin size, with the 
objective to create bins as independent from each other 
as possible. Large Cherenkov detectors use sparse sam- 
pling with sizeable distances between sensors, and each 
sensor observes light over a limited track segment. The 
length of the segment depends on the optical proper- 
ties of the medium, with different properties for sea wa- 
ter, ice, and liquid scintillators that vary with the wave- 
length of light. 

In water, absorption effects are much stronger than 
scattering. Thus the most appropriate measure of op- 
tical properties would be the absorption length, which 
peaks at wavelengths ~470 nm. The absorption length 
is ~22 m in lake water 13211 and ~60 m in seawater 13311 . 

In deep glacial ice, however, scattering effects are 
stronger than absorption, and the appropriate measure 
is the light propagation length, which peaks ~400 nm. 
This length depends on both the effective scattering 
length A es and the absorption length A a in the medium, 
using the relation 



(4) 



In IceCube, for example, the average absorption co- 
efficient is ~0.008 m , or an absorption length of 
125 m, while the average effective scattering coefficient 
is ~0.040 irT 1 , or an effective scattering length of 25 m 
13. The average propagation length is thus ~32 m, 



which increases to ~45 m for the very clear ice in the 
bottom of the detector. 

For liquid scintillator detectors, the absorption length 
is the appropriate measure of optical properties, which 
is typically 10-20 m fl^l. 

The detector geometry may also impose constraints 
on the truncated mean's binning method. In ANTARES, 
for example, the string spacing is 70 m, which is 
more than the reported absorption length. In Hyper-K, 
with its large photocathode coverage, the water proper- 
ties may be the determining factor for the optimal bin 
length. Hyper-K's 250-m tank length is far longer than 
the light attenuation length, and the detector geometry 
will naturally allow for track segmentation. 

Most of the vertical strings of optical modules in Ice- 
Cube are deployed on a 125-m grid, making the detector 
roughly hexagonal in shape. For near-horizontal tracks, 
this string separation sets the minimum sampling dis- 
tance, and hence a minimum bin length. In this study, 
since the IceCube array is not perfectly regular and to 
avoid having or 2 strings in a bin for horizontal tracks, 
a 120-m bin length was used, spreading the energy more 
evenly across the bins. Due to the size of the detector, 
events could have up to 15 bins. 

Figure [4] shows a scatter plot of simulated versus 
truncated dE/dx for the same events as in Fig. [TJ but 
now using the truncated dE/dx with 40% of the high- 
est bins discarded, and including the optimizations dis- 
cussed in Section [3] For this analysis, a minimum of 
3 bins was required for the event to qualify, since a 
40% truncation would result in bins omitted for a 2- 
bin event. The distribution is much narrower than the 
conventional untruncated dE/dx method but retains the 
same S-shape, for the same reasons. 

The spread in dE/dx values is much narrower for the 
truncated method than for the conventional method. A 
comparison of the distributions of dE/dx values for the 
same narrow energy range as in Fig.|2]is shown in Fig. [5] 
For the truncated method, the Gaussian cr is 0. 18, com- 
pared to the untruncated value of 0.29, which constitutes 
a significant improvement, and the skewness is much re- 
duced. 

Figure [6] shows the corresponding energy resolution 
using the 3-equation fit from Fig.|4]to calculate the en- 
ergy. Compared to Fig. [3] the high-energy tail is much 
smaller in the truncated plot, with only 1.3% of the 
events reconstructed with more than 5 times their ac- 
tual energy. The Gaussian resolution has improved to 
0.22 in log 10 (£' /i ), which is a factor of 1.6 as compared 
to 2.0 for the conventional untruncated method (with the 
3-equation fit). 



6 



10r 





£ 3000 : 



1 2 3 4 5 

log (Truncated dE/dx [GeV/m]) 



Figure 4: (color online) Simulated vs. truncated dE/dx at the cen- 
ter of the detector, for the optimized 40% truncation. The distribution 
is much narrower than in Fig.fJJ The vertical dashed lines identify the 
ranges for the three fit equations (represented by the solid black line). 
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Truncated dE/dx 
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Figure 5: (color online) Truncated muon dE/dx vs. conventional (un- 
truncated) muon dE/dx, for simulated energies ~10 TeV (log^i?,,) 
= 4.0 ± 0.05). For the truncated method, the Gaussian cr = 0.18 ± 
0.01, and RMS = 0.22, compared to the untruncated method values of 
cr = 0.29 and RMS = 0.44. The effect of the large stochastic events 
has been greatly diminished, thus reducing the spread. 



3. Algorithm Optimization and Performance 

We studied several variations of the algorithm to op- 
timize performance. The parameters included the frac- 
tion of bins discarded in the truncated mean (Section 
13. Il l and limits on the distance between the sensor and 
the reconstructed track (Section [3~!2l i. We also evaluated 
the inclusion or exclusion of the unhit sensors in the cal- 
culation (Section [3~3l ). and we explored a variation in the 
truncated method that treated each sensor as its own bin 
(the "DOMs method") (Sectionl3~4l). Lastly, we looked 
at the median versus the mean for both the binning and 
DOMs methods (Sectionl33i 



1000 




1 2 
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Figure 6: (color online) Truncated muon energy resolution, using the 
fit curve from Fig. f4] for muons with simulated energies between 
1 TeV and 1 EeV. The number of entries to the right of the Gaus- 
sian curve (black line) is much smaller with the truncated method, as 
compared to Fig. [3] The energy resolution improved to cr = 0.22 in 
log )0 (£^) (a factor of 1.6 in energy), with the mean of the Gaussian 
curve slightly offset by -0.03, and an RMS value of 0.25. 



3.1. Truncation Percentage 

The optimal truncation was determined by evaluat- 
ing truncation percentages between and 80% of the 
bins. If the percentage of cut bins resulted in a non- 
integer value, the fraction was rounded down (for exam- 
ple, 40% of 4 bins would result in 1 bin cut). Figure [7] 
shows the resulting energy resolutions for different trun- 
cating percentages (up to 60% shown). In this analysis, 
40% was the optimal cut. 
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Figure 7: (color online) Energy resolution for truncation percentages 
ranging from 0% to 60%, for E^ > 1 TeV. The optimal truncation 
percentage was 40%. Gaussian cr were 0.27, 0.23, 0.22, and 0.23 
for 0, 20, 40, and 60%, respectively. The optimization techniques 
described in Sections l3.2l and r3~3l have already been applied. 
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3.2. Distance of optical modules from track 

Energy measurement by dE/dx relies on an accurate 
measurement of the distance between the tracks and the 
optical modules. However, the reconstructed track po- 
sitions (transverse to the direction of motion) can be off 
by several meters. Also, data from optical modules that 
are very far from the track are less useful because of the 
increased chance that the sensor hits are noise. 

To avoid these issues, we evaluated the performance 
of the method when requiring a minimum and maxi- 
mum distance between the track and the sensors. The 
resolutions for various distances are shown in Table Q] 
for horizontal muons with simulated energy > 10 TeV. 
The optimal energy resolution used optical modules that 
were located between 10 and 80 m from the track. In an 
optically noisy environment, such as seawater, a smaller 
maximum distance likely would be optimal. 





Cylinder distance (meters) 


Cuts 


None 


0-100 


10-100 


0-80 


10-80 


0-60 


10-60 


0% 


0.362 


0.345 


0.331 


0.342 


0.328 


0.329 


0.308 


10% 


0.324 


0.316 


0.282 


0.307 


0.291 


0.319 


0.298 


20% 


0.276 


0.273 


0.222 


0.244 


0.216 


0.242 


0.221 


30% 


0.255 


0.231 


0.211 


0.228 


0.206 


0.229 


0.204 


40% 


0.248 


0.227 


0.210 


0.229 


0.198 


0.225 


0.204 


50% 


0.250 


0.234 


0.215 


0.229 


0.200 


0.227 


0.202 



Table 1: Energy resolution cr in log| (£', 1 ) for varying truncated per- 
centage cuts and cylinder distances (in meters) from track, for muons 
with simulated energy > 10 TeV. The optimal resolution occurred with 
40% truncation and a cylinder of 10 to 80 m. When the lower-energy 
muons are included (E„ > 1 TeV), the energy resolution is slightly 
degraded (<x = 0.22). 



3.3. Optical modules without observed photoelectrons 

In IceCube, most energy reconstruction methods 
(such as the truncated mean method) use a Gaussian 
error model, which is valid at high energy but not at 
low energy where the number of detected photons is 
bounded at the low end by zero but is not bounded 
above. For optical modules that detect many photons, 
this is mostly irrelevant because Poisson errors become 
symmetrical. Therefore, the inclusion of many low- 
amplitude sensors, such as those that did not observe 
any photons (the "unhit" sensors), could bias the likeli- 
hood calculation. 

In addition, isolated hits (optical modules without 
hits in their nearest or next-to-nearest neighbors) are 
likely to be noise. For these reason, these energy-loss 
calculations only used the data from optical modules 
where the nearest or next-to-nearest neighbors were hit 
within + 1 /is |35|, otherwise known as coincident-hit 
modules. 



In calculating the number of expected photoelectrons, 
there is no distinction between hit and unhit optical 
modules, and Photonics will return non-zero amplitudes 
for all optical modules. To reduce the effect of this mis- 
match, we studied the effect of counting the expected 
photoelectrons from only the optical modules with coin- 
cident hits, rather than using all of the optical modules. 

Figure[8]shows the results from this alternate method- 
ology combined with the truncated mean. The energy 
resolution for this "hit" technique was slightly better 
at the higher energies. However, at lower energies, the 
curve is much steeper than in the previous "all-sensor" 
truncated curve, and the energy resolution was signif- 
icantly worse when using only the hit sensors in the 
dE/dx calculation. Thus, this "hit approach" was not 
used in the final truncated method. Instead, the data 
from the coincidence-hits (within the cylinder) were 
used for the actual photoelectron counts, and all the sen- 
sors (hit and unhit) in the cylinder were used for the cal- 
culation of expected photoelectrons. 




log i0 (Truncated dE/dx [GeV/m]) 



Figure 8: (color online) Simulated vs. truncated dE/dx, using only 
the "hit" optical modules within the cylinder. The vertical dashed lines 
identify the ranges for the three fit equations (represented by the solid 
black line). At low dE/dx values, the curve drops off more steeply 
than in Fig. [4] which adversely affects the energy resolution in this 
range. 



3.4. DOMs method 

An additional option was explored that diverged from 
the binning method, which was nicknamed the "DOMs 
method" after the digital optical modules in IceCube. 
In this option, instead of using 120 m bins, each op- 
tical module became a "bin" and was given its own 
dE/dx value by dividing the observed number of pho- 
toelectrons by the expected number of photoelectrons. 
The advantages of this method are (1) the track can be 
shorter (no minimum number of bins), (2) additional 
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data will be used in the energy determination, and (3) 
the algorithm is slightly faster. The main disadvantage 
is that a large stochastic event would be recorded by 
many sensors, and thus there would be an undesirable 
correlation between sensors. 

The technique was optimized (using similar proce- 
dures to those outlined above) to require a minimum of 
8 optical modules within to 60 m of the reconstructed 
track for the event to qualify (versus 10 to 80 m for the 
bins method), and the optimal truncation was to omit the 
highest 50% of the optical modules (versus 40% of the 
bins). Then the dE/dx values for the remaining optical 
modules were averaged. This yielded an energy resolu- 
tion and RMS improvement that was similar to that of 
the binning method (cr = 0.22 and RMS = 0.25). 




1 2 3 4 5 

log (DOMs Method dE/dx [GeV/m]) 



Figure 9: (color online) Simulated vs. truncated dE/dx (DOMs 
method) at the center of the detector. The vertical dashed lines iden- 
tify the ranges for the three fit equations (represented by the solid 
black line). The event distribution is much narrower than in Fig.fT]but 
similar to the binning method in Fig. [4] 

Figures[9]and[T0]show the results of this method. One 
or the other method may prove more useful in the var- 
ious detector designs and is an option that can be ex- 
plored. In this analysis, about 4% of the events would 
qualify for only one of the methods, while the remainder 
would qualify for both. Including the two approaches 
maximizes the capability of the truncated mean method. 

3.5. Median Method 

One final option was to change from the truncated 
mean to the median of the bins. In this method, the 
truncation occurred at both the upper and lower ends 
of the ordered list of bins (or DOMs, in the case of 
the DOMs method), such that only one bin (DOM) re- 
mained. The dE/dx value for this bin (DOM) became 
the value for the event. For both the binning method and 
the DOMs method, which had separate optimizations, 
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Figure 10: (color online) Truncated muon energy resolution for the 
DOMs method, for muons with simulated energies between 1 TeV 
and 1 EeV, using the fit curve from Fig. [9] The number of events to 
the right of the Gaussian curve is much smaller than in Fig. [5] The 
resolution is very similar to that of the optimized binning method as 
shown in Fig.[6](cr = 0.22 and RMS = 0.25). 



the energy resolution and RMS were nearly unchanged 
for all muon energies when compared to the earlier trun- 
cated mean methods. 



4. Results and Conclusions 

There is no apparent correlation between the energy 
calculation and the zenith angle (see Fig.fTTTi. The distri- 
butions of the energy residuals are very similar in shape 
and width. There is also no apparent correlation be- 
tween the energy calculation and azimuth angle, optical 
properties of the medium, or other parameters. Thus the 
truncation method is universally applicable to all par- 
ticle track zenith and azimuth angles within the detec- 
tor, with the proper bin size. The energy resolution im- 
proves with an increasing number of bins as expected, 
as shown in Fig. [12] but levels off at 0.18 in log^E^). 

Figure [13]compares the RMS and Gaussian cr values 
of truncated and untruncated dE/dx for various muon 
energies, showing the improvement from using the trun- 
cated method. 

Figures [141 and [TBI contrast the untruncated and trun- 
cated methods, respectively, for actual versus calculated 
muon energy. The improvement in energy resolution is 
fairly uniform over the energy range. 

Figure [16] compares the input muon energy spectrum 
to the spectrum of calculated energies from the trun- 
cated method. The agreement is quite good, with a 
slight shift in energy below 1 TeV due to the off-peak 
fit equation plus the spread in dE/dx at low energies. 
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Figure 11: (color online) Energy resolution (in log 10 (£ /J )) for zenith 
angle bins of 30 degrees for from 1 TeV to 100 TeV, using the 
truncated bins method. There is no visible correlation between zenith 
angle and energy resolution, for the energy range from 1 TeV to 1 EeV 
using an E~ spectrum. 



Figure 13: (color online) Comparison of the RMS and Gaussian <x 
of untruncated and truncated dE/dx values vs. log 10 (£ , , 1 ), for simu- 
lated energies between 1 TeV and 1 EeV. The energy resolution for 
the truncated method is much improved over the untruncated method. 




10 12 14 
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Figure 12: (color online) Correlation between the number of bins and 
the average energy resolution (in log 10 (£', 1 )) of the muon events, using 
the optimized truncated mean method, for E^ > 10 TeV. 



However, there are no glaring discontinuities in the cal- 
culated spectrum, which indicates that the truncated 
method gives back the original spectrum quite well. 

By using the truncated mean method, the energy res- 
olution is significantly improved. The best truncated 
method incorporated the following criteria: (1) only in- 
clude the photoelectrons from optical modules within 
10 to 80 m of the track, (2) truncate the highest 40% of 
the bins, (3) use both hit and unhit optical modules in 
the calculation, and (4) sum the remaining photoelec- 
trons separately (observed and expected) to determine 
the new truncated dE/dx value. With these optimiza- 
tions, the energy resolution was improved from 0.29 
in log 10 (£^) to 0.22, for the energy range of 1 TeV to 
1 EeV. This is a 26% improvement in the overall energy 
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Figure 14: (color online) Actual muon energy versus calculated muon 
energy for the untruncated method with the 3-equation fit. The black 
line is a perfect 1 : 1 correspondence. 



resolution, with better resolutions above 10 TeV. The 
technique is applicable to any detector that uses dE/dx 
as the primary means of the particle's energy determi- 
nation. 
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Figure 15: (color online) Actual muon energy versus calculated muon 
energy for the truncated method (40% cuts with optimizations). The 
black line is a perfect 1 : 1 correspondence. 
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